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Abstract 

We use the transfer matrix formahsm to derive non-perturbative sum rules in Wilson's 
lattice QCD with A^f flavours of quarks. The discretization errors on these identities are 
treated in detail. As an application, it is shown how the sum rules can be exploited to 
give improved estimates of the continuum spectrum and static potential. 



1 Introduction 



Although relatively old, the subject of lattice sum rules ^ in d = 4 SU(A^) gauge 
theories has received only sporadic attention since it was initiated by C. Michael In 
essence, these identities relate the derivative with respect to a bare lattice parameter of 
a physical quantity to a higher order correlation function. A prototype sum rule states 
that the /3-dependence of an energy level is determined by the expectation value of the 
Lagrangian on that state. 

In Monte-Carlo simulations, sum rules were first applied as a cross-check in form 
factor calculations at zero- momentum transfer [1], and in measurements of the flux- 
tube profile in the static QQ system 0. Later they were used to determine lattice 
beta-functions non-perturbatively jHllZI- On the theory side, they have also led to some 
insight into the relative contribution to the mass of a glueball of the component Tqo of 
the traceless energy-momentum tensor and of the trace anomaly jSj. 

Here we shall generalize the sum rules to the case where A^f flavours of Wilson quarks 
are present. We do so by invoking the transfer matrix, in our view the most elegant 
approach. It also makes it clear that the identities derived do not depend on boundary 
conditions or on which phase the system is in — be it a finite-temperature phase, a 
parity-broken phase or other. Let us recall the Wilson action, which has two parameters 

5g = PY^Up, □p = ^Re'IV{l-C/(p)}, (1) 

p 

Si = Yipix)ipix) - nYlfj. V'Wf^AiC^;)!! - 7m)V'(^+«A) + V'(^+»A)f^MW(l + 7Ai)V'W (2) 

X 

where in Eq. ^the sum extends over all unoriented plaquettes. 

It is well-known that cutoff effects are of 0(a) in Wilson's formulation. In recent 
times, most large-scale numerical calculations have made use of the 0(a) improvement 
program to reduce these effects (with some notable exceptions jlH I12j). We present 
an idea to exploit lattice sum rules in order to extract additional information on the 
lattice spacing dependence of the spectrum, and further how its lattice artefacts can be 
reduced. 

The idea is also applicable to the energy levels in the presence of a QQ pair, where 
one is typically interested in their derivatives with respect to the static quark separation. 
The static force, and the effective central charge in the pure gauge case, are affected by 
discretization errors which can be large compared to the statistical precision one is able 
to achieve. 

In section 2 we rederive some of the known results in the pure gauge theory, where 
particular attention is paid to discretization errors. In section 3 the sum rules are 
derived for full QCD. Concrete applications are proposed in both sections. We make 
some concluding remarks in section 4. 



2 Transfer matrix and sum rules in the pure gauge theory 

We start with a brief reminder on the transfer formalism and how sum rules are derived 
in the pure gauge theory. Some familiarity with the lattice regularization is assumed. 
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For dimensionful quantities in lattice units we use the notation L = La etc. 

If V is the set of spatial link variables in a time-slice, let ^\y] £ Tic be a square- 
integrable wave function with respect to the Haar measure of the gauge group. The 
transfer matrix acts as follows I14| 113]: 

i-r^nV] = [U^LidV^i^] Kf,[V,V']^Vl (3) 
with kernel Kp[V,V'] = / JJ dVF(x) exp (-/3A5[y, T^, F']) , (4) 

X 

AS[v,w,v'] = Y.{^kOok[v,w,v'] + hEk<iiOki[v] + aki[v'])]{5) 

X 

nok[V,W,V'] = ^ReTr{l-14'(x)VF(x+aA;)yfe(x)-ipF(x)-i}. (6) 

In words, AS is the restriction of the action to the region of space-time between xq and 
xq + a, with a weight of 1/2 given to the terms living on the boundary. When the full 
gauge system is defined with periodic boundary conditions in the time direction, the 
partition function is given by Z(P) = TrlT^^"}. 

We define the magnetic-plaquette operator Clki to act as in Eq. El with Kjj replaced 

by 

Kki [V, V] =6{V-V') Oki (0, x) [V] . (7) 
The kernel corresponding to the electric plaquette operator Dofc is defined as 

Kok[V, V']= / J] dl^(x) exp {-pAS[V, W, V']) DokiV, W, V']. (8) 

X 

Then, if T|$) = A|$), 

-dpxip) = -d^{<P\Tf,\<P) = -{^\dpT0\<^) = x{p){^\Zk<Pkm + mZk^okm- 

(9) 

Since |$) can be any eigenstate of the transfer matrix, we can choose it to be successively 
the vacuum and a glueball state. The spatial volume is assumed to be large enough 
for the finite-volume effects on the mass gap to be negligible. Subtracting the two 
corresponding rules gives 

^ = Em(^Ph'J>) - i^Pkin) + Efce*'(^Pofc|^> - e^^{n\\Jom. (10) 

Here M = — log A/Aq is the mass of the glueball in lattice units. One might be surprised 
to see the vacuum energy appear in this expression, but as we shall see shortly, computing 
the expectation value of an electric plaquette in the usual ensemble yields precisely 

Following T!, we can use the fact that d/ da{M / a) = in the continuum limit to 
conclude 

(11) 
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To translate the matrix elements appearing on the right-hand side into Euclidean 
correlators, for $ the lightest state in its symmetry channel one may choose a zero- 
momentum linear combination 4'{xo) of magnetic loops in that channel, and compute 
the correlator 



(^(xo + a) i ( Efc</ Pkiixo/2) + Pki{xo/2 + a)] 0(0)) 



^{<^>\a,im, (12) 



(</.(xo + a) m) 

where (j){xo) = (j){xo) — {(/)) (take xo/a even). And similarly 

(0(xo + l) EkPok{xo/2) 0(0)) .0^^ e'^''^('^S^(<^\aM (13) 
(0(xo + a) 0(0)) ' °" ^ ' 

2.1 Improved estimators of the continuum spectrum 

We use the shorthand notation = for Eq. ^Jfrom now on. In view of taking 
the continuum limit, we use a reference length £q which sets the scale of the theory. 
Typical examples are o""^/^, where a is the string tension and ro, the Sommer reference 
scale (l£j. We assume that for every £q = io/a can be determined, so that P{io), 
which is assumed to be monotonic in the range of interest, and 

d(3 _ dP{£o) ^^^^ 



d\oga dlogio 

are well-defined. For the case £q = ro, ioiP) is known with one percent precision or 
better jJH]. Thus on a line of constant physics, a function of f3 can just as well be 
regarded as a function of io, and vice versa. 

The dimensionless quantity z = IqM has a continuum limit, and we can express its 
derivative with respect to the lattice spacing exactly: 

t = M--^nt (15) 

dio dloga 

This information on the slope can be included in the continuum extrapolation if data at 
several lattice spacings is available. Even if not, it can be used to provide an improved 
estimate of the continuum limit of z if we assume a particular form of the discretiza- 
tion errors. In practice z is usually observed to approach the continuum with O(a^) 
corrections within statistical errors; such corrections can be removed according to: 



impr 



-(4) + k^o^ = Ui,) - at (16) 

d£n 2 2 dloga 



In what sense is this an improved estimate of the continuum limit? Suppose the true 
analytic form of z is [22 IHj 

1 ^ 

ziio) = W + )i E log"{l/4} + 0(1/^"^). (17) 

^0 n=0 
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Then 

1 

^impr = Zcont " + l°g"{V^o} + 0(l/£^). (18) 

^^0 n=0 

Since the series of logarithms is asymptotic, if = 1 yields the best accuracy at the 
lattice spacing one is working at, then |co| ^ |ci log{l/£o}| ^ |ci|, so that zimpr has 
reduced discretization errors, although it is not a full O(a^) improvement. Note that it 
is essential in (|16|) to use the quantity £q that appears elsewhere in the formula to define 
A, as in Eq. d- 

2.2 States depending on an external length scale 

We now consider an eigenstate of a transfer matrix which, in addition of depending on /3, 
also depends on a physical length scale £, and the latter can take only integer multiples 
of some quantum As: 

s = n As, n e N. (19) 

Examples are zero-momentum states in finite volume (where s = L, As = a), but also 
large- volume states with non- vanishing momentum (s = p, As = 2tt/L). The sum rule 
then reads 



—in,p) = a-''''^a-^. (20) 

We shall make use of the standard notation for discrete differences 



dnfin) = /(n + l)-/(n), d*Jin) = f{n) - f{n - 1) (21) 
dnf{n) = i(/(n+l)-/(n-l)), A„/(n) = /(n + 1) - 2/(n) + /(n - 1). 

Let us consider again the quantity z = (.qE, which has a continuum limit. For 
definiteness take the case of a zero-momentum state in finite volume, so that L = na. 
For a given couple (n, /3), we choose an auxiliary n' close to but different from n; there is 
a /?' such that n'/n = £o{(3')/£o{f3); that is, the box sizes are matched in physical units. 
Then we can write 

ioE{n,P) = £o{(3')E{n',p')[l + 0(0^)] (22) 



^0^ 



Ein'^f^) + l&-a^°s{n/n'}n^i(3) + ^ log'' {n/n'}f {n' , f3) [1 + 0(0^)] 



where /(n',/3) = dp E{n' , P) + 8^^ E{n' This second order term is 

necessary to eliminate 0(a) cutoff effects. We now choose n' = n — 1 and then shift 
n — > n + 1, obtaining Eq. (a); secondly we choose n' = n + 1 and then shift n — > n — 1, 
obtaining Eq. (b). Note that the function f{n,P) is evaluated at the same arguments 
in (a) and (b). Therefore we can eliminate that term by taking the linear combination 
log2(l - 1/n) • (a) - log2(l + 1/n) • (6). This results in: 

A an [1 + o(a2)] = Bin) + ndnEin) (23) 

where we dropped the (3 dependence and used the freedom to trade one discretization 
scheme for another that is equivalent to O(a^). This is one of the Michael-Rothe sum 
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rules ^ 12] . We have however kept track of the discretization errors carefully so as not to 
introduce 0(a) effects. This sum rule allows one to extract the derivative of the torelon 
energy with respect to its length without having to perform an independent simulation. 
In particular the effective central charge can be computed via 

Ceff(L = na,£o) = h^a^l - ^E{n). (24) 
2.3 The static potential case 

The sum rule (|20|) also holds if static charges are inserted at points and x, x = 
rk. Indeed the kernel of the transfer matrix then projects onto states ^ai3... which, 
under a gauge transformation U — >■ U^, transform with A(0) and A(x) according to 
the representation of these charges. For a fundamental-antifundamental pair, we have 
explicitly: 

^a(3[U^] = {HO)a^T A(x)^5 ^^s[U]. (25) 

The transfer matrix kernel is 

^n'^In = / n W W, V']) {W{OUr W{^)^s (26) 

and the composition rule is 

(T^*0n)"^^' = (Tn*«n)""^^ (Tn.^n)"^^''. (27) 
Thus the partition function in the presence of the static charges is 

^N*^N = ]^ E {(Tfc^N)"""^ (28) 

a, 7 

In particular, the Polyakov loop correlator evaluates to (-P(O)* -P(x)) = Zj^^^s^/Z. 

The energies at two different separations r must be subtracted in order to remove the 
divergent self-energy of the static quarks. The quantity £o(/3)5„£'(n, /?) has a continuum 
limit. We define n = n + ^ -|-0(a^) where the O(a^) need not be specified presently, and 
choose two integers n,n' ^ 1, \n — n'\ = 0(1) and a particular value of /3. Now there 
exists an auxiliary /?' such that 

na{P) =^a{/3'), i.e. nio{p') =^£o{(3). 

We can write £l{p)dnE{n, /3) = £l{(3')dn'E{n', (3')[l + 0(0^)] and then Taylor-expand in 
/? to obtain in the same way as in the previous section 

nAnE{n) + 25„^(n) = ^^^^tl [1 + 0(0^)]. (29) 

We use the notation D"^]; = Dq^^ — Dj^- Again a symmetric finite-difference scheme is 
necessary and sufficient to remove the 0(a) discretization errors. 
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2.3.1 An improved estimate of rQF(r) 

We now show how one can also use the lattice sum rule to reduce cutoff effects on the 
static force. Discretization errors are by far the dominant source of uncertainty on this 
quantity P^ . 

to = ro and z{^, fo) = rlF{r = ^tq, ro). (30) 
A direct measurement of the force is 

F{na,ro) = dnE{n,ro). (31) 

This define the static force at a discrete set of points. To define it for all distances, an 
interpolation formula between neighbouring points must be used. Which formula is a 
matter of choice (different definitions of the force at finite lattice spacing will then differ 
by 0(a2)). 

For illustration we choose a linear interpolation between the two nearest direct mea- 
surements of the static force: 



^2 

z(^,fo)- ° 



n — n — 1 



{(ro -n- l)dnE{n, ro) + (n - CWn^in, ro) 



(32) 



With this precise definition and Eq. [201 dz/dro can be evaluated at fixed ^ exactly. An 
improved estimate of the continuum z reads 



2impr(C, ^o) = z{(, ro) + h^o-QU^, (33) 



which is improved in the same sense as the glueball mass in section I2TT1 Explicitly, using 
Eq. EOl we get 

1^2 

^impr(e,ro) = 2z{(,ro) + ^^J= X (34) 

n — n — 1 

(roAnEin, ro) - [m - 7r^)n^+' + (n - (ro)a^-i 

It is important to realize that the improvement is not merely a higher order difference 
scheme for the static force, rather it contains non-perturbative information about its 
lattice spacing dependence. Formula (|35j) may be numerically useful since all terms 
on the right-hand side can be evaluated in the same simulation. Note that the idea is 
applicable to more complicated interpolation schemes for any particular definition of n, 
and also to the effective central charge Ces{r) = —^r^cPV/dr'^. 



2.4 Anisotropic couplings 

Sum rules in the pure gauge theory derived by varying the lattice spacings in the four 
space-time directions independently around the isotropic point can be found in We 
have little to add to this subject, except to say that the sum rules involving an external 
scale, such as a lattice size dependence, must be expressed with a symmetric finite- 
difference scheme to avoid 0(a) discretization errors. Here is a nice application: for a 
torelon in the direction 1, the energy density in the transverse plane is directly sensitive 
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to quantum string corrections, as already noted in and an effective central charge 
can be defined via 



n 



Ces{L = na,£o) = t; {U - S) 



□02 + □ 



03 



2n. 



23 



(35) 



The definition of the anisotropic derivatives U and S will be given in section 13. 2t we 
note that {U — S)/2 = (3 — N{ca — Cr), where the coefficients Ca, Cr were determined 
in (see also Ref. therein) for a few /3 values. The effective central charge extracted 
in this way will differ by O(a^) terms from Eg. 1241 



3 Lattice sum rules with Wilson fermions 

We start by recalling some essential facts about Wilson fermions and their transfer 
matrix. We shall use x, y for three-component spatial vectors, indices a,/3, . . . for color 
and p,a, . . . for spinor indices. A spinor component is written as V'(t(x, a); for simplicity 
we consider the one-flavour theory for the moment, but the extension to several flavours 
is trivial, as we shall see. We use a set of Euclidean Dirac matrices, {7/i,7t/} = 25 f^^,, 
as well as the projectors P± = ^(1 it 70). Whenever there is a risk of confusion, we 
use a hat to distinguish a quantum mechanical operator from a Euclidean c-number or 
Grassmann variable. The full Hilbert space Tl is now the tensor product of the gauge 
Hilbert space Tic and a fermionic Hilbert space fip- 

An important object is the color-covariant, nearest-neighbour transport operator for 
fields transforming in the fundamental representation: 

D^{^,a;y,p) = Uk{x)ap 6^^',y ± Ul{y)ap 6y^f^^^, A; = 1,2,3. (36) 

Two spatial finite-difference operators will be used: 

B = 1-kY,D+ (37) 

k 

C = \Y.^k®lk. (38) 

k 

B is hermitian and strictly positive for < k < 1/6, while C is anti-hermitian , for any 
gauge field configuration [H] . 

The fermionic Hilbert space is the Fock space built from a collection of x operators, 
which enjoy canonical anticommutation relations {{Xai^iOi), x^'(y)/3)} = '^xy ^aa' ^ap 
etc.). One further defines Jl] 

= = £ {B-^l^)\ (39) 

that turn out to be the operators associated with the Grassmann variables (-0, ■(/')■ 

The transfer matrix in the temporal gauge ^0 = was obtained in terms of the x 
operators in |14lll5j . The kernel of the transfer matrix in the path integral form is given 
in [201, Eq. 4.26: 

K[V,i^,^P,V',i^',^'] = det{5fi'}i/2 /■-Q^^(^) g-AS[y,w^,y'] ^40^ 

exp{-2(ViP+K;CV' - VSP+kH^" V - i^'WP-K^l) + tp' P^kC'^j')}, 
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where B', C' are functionals of V'. The scalar product of two functionals of ('0, ip) is 
taken with respect to Grassman integration with a measure e~^^^/detB. The deter- 
minant arises because of the change of basis (|39|) . 

We can now derive a sum rule in the same fashion as in the pure gauge case. Let $ 
be a normalized eigenstate of the transfer matrix with eigenvalue A. We find^ 

d^X = 5«($|T(/3,k)|$) = 2 j DVDijD^P det{B}-^/^ e"'^^^ ^>[y,Vi,V]* 

I DV'Dij'Di^' det{S'}-i/2 e-^'B'V'' $[^',^^',^'1 

ex.p{-2{tpP+KC^p - i^P+KW^i)' - tp'WP-Kip + tp' P^kC'^j')} 

DVDijDij -B)- 2C)V) (41) 

In the last line, A multiplies the expectation value on the state $ of an equal-time 
operator ^1^, while the rest of the expression is the expectation value on the state $ of 
an integral operator Af. Following the discussion of the pure gauge case, we thus have 

-d,E{p,K) = - {n\As\n) + e^{<^\At\<^)-e^"{n\At\n) = h\t (42) 

To evaluate this matrix element in a Monte-Carlo simulation, one would use an inter- 
polating field if for the state If it is the lightest in its symmetry channel, 

($|i,|«>)= lim fsixo), e^mAt\^)= lim ft{xo) (43) 

with (for xo/a even) 

f ^„ \ - {'?(0)^{As{xo/2)+A4xo/2+a))^(xo+a)) ^ \ _ {(p{0)At(xo /2)(p{xQ+a)) 

Js{XO) - {^{0)<fi(xo+a)) ' - {<p{0)<fi(xo+a)) y^^> 

and 

Asixo) = ^Vi(x)C/fc(x)(l-7fc)V'(x + aA:) + V5(x + aA:)C/fc(x)~i(l+7fc)V'(x)(45) 

x,A; 

At{xo) = ^V^(x)f/o(x)(l-7o)V(x + aO) + V5(x + aO)(l+7o)[/o(x)-V(a;).(46) 

X 

These are the hopping terms of the Wilson-Dirac action; they are the terms one would 
obtain by naively differentiating the Boltzmann factor. It is also clear now that equa- 
tions (|^T|) and (|l5lll6|) hold for A^'f degenerate flavours, provided an implicit summation 
over flavours is understood in ()451 I46() . In practice, once the fermion fields are inte- 
grated out, the three-point functions involve all-to-all propagators within one time-slice. 
A stochastic estimator is then required, as used in previous thermodynamics applica- 
tions |24j . 

""^It may be useful to note that if T{x,y) = T*{y,x), \ — f w{x)dxw{y)dy(f>* {x)T{x,y)4>{y) 
and ^ (j)* [x)4){x)w{x)dx — 1, where T, (j) and w all depend on a parameter k, then 9^ A = 
J 'w{x)dxw{y)dy(l>*{x){d^T{x,y))(t){y) + \ J{d^w{x))dx\(t){x)\^. 
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3.1 Applications 

We now present some applications of the sum rule just established. We focus on the case 
of A^f degenerate flavours. Let us assume for now the spatial volume to be large enough 
for the state $ to suffer negligible finite volume effects (for a stable one-particle state, 
they are exponentially small). Given the existence of the continuum limit, we have, up 
to cutoff effects. 







d_E 

da a 



(/3(a), K(a))] 



E = dpE 



d(3 



dlog a 



+ d^E 



dn 



dlog a 



(47) 



But actually, k = n[l3{a)): on a line of 'constant physics' n must be tuned as a function 
of 13. Hence 



E 



dp 



dK (J, 



[1 + 0(g)]. (48) 



dlog a 

where we have used Eq. EH and Eq. |12l Recall that j^^^ 
continuum limit (6o = 3(4^)^ (llA'" — 2A'f)). A particularly interesting case arises for 
iVf > 2 in the chiral limit, k = Kc- The 'pion' mass then vanishes and hence 



4:Nbo is universal in the 



^ = S|-[l + 0(a) 
dp h\f^ 



(49) 



Note that at no stage in the derivation of the sum rules did we need to make an as- 
sumption on the extent of the time direction - only the spatial lattice size was assumed 
to be essentially infinite. Hence this equation also holds in the Schrodinger functional - 
where one can actually simulate at Kc- Relation (|49)) can now be inserted into the sum 
rule for the only other stable particle, the nucleon: 



nucl 



dp 



□ 



nucl 



1, 1 nucl 



[1 + 0(g)]. 



(50) 



dlog a 

As in the pure gauge theory, we can give an improved estimate of the continuum 
z = IqM, assuming now that the leading corrections to the spectrum are 0(a). Here 
■^o(/?) f^ciP)) is the quantity evaluated in the chiral limit; concrete examples are ro|m^=0) 
A^ or Lmax) defined by 5^(l/imax) = x, where x is a particular numerical value and 
(p' the renormalized Schrodinger functional coupling j23j . We get 



2impr(4) = Z{k) + dz/dlogk = 2z{io) - 



dp / $ dK ^ 



(51) 



We expect the 0(a) cutoff effects to be substantially reduced on this estimator. Note that 
for this formula to yield an improvement, must be defined through the quantity 

io, as in Eq. (|TH) . 

If the volume dependence of the energy level cannot be neglected. Eg. 1481 becomes 



Ei{Lk}) + ^Lkd^L, 

k 



E 



dp 
dlog a 



^ dK , 



:i + 0(a)]. 



(52) 
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The derivation is identical to the one for the corresponding relation in the pure gauge 
case. The finite-difference scheme used here in the derivative is formally irrelevant, since 
the sum rule holds with 0(a) corrections anyhow; however we still expect the symmetric 
scheme to yield somewhat smaller discretization errors. This formula may be useful to 
study the volume dependence of the pion mass, and also of the state which is a mixture 
of a two-pion state and the p resonance. The latter volume dependence allows one to 
extract the scattering lengths of pions in the elastic regime and the p width [22j . 

3.2 Anisotropic lattice 

We consider the Wilson theory with A'^f degenerate flavours. If one assigns an indepen- 
dent hopping parameter in the space directions from the time direction, then one can 
derive the equations 

e^{^At\^) - e^^{n\At\n) = ht\% (53) 
($|i,|$)-(J7|i,|f7)^h,|H (54) 

To exploit these relations in a useful way, we must understand how the theory with 
parameters (/?s, /9f , Kg, i^t) is renormalized. As independent 'physical' variables, the most 
intuitive choice is (^O)'^)^;)) where ^ is the anisotropy ag/at and x essentially sets the 
quark masses (e.g. x = mpg/my). The continuum limit is then taken with £g ^ — > at 
fixed and x. One can alternatively choose the set (of, independent variables. 

Since there are four lattice parameters, it is only on a hypersurface in that parameter 
space that the lattice theory describes QCD at all. The continuum limit corresponds to 
a particular curve on that hypersurface where ^ and x are constant. 

Imagine momentarily using four independent lattice spacings, so that (ao, ai, Q^3) x) 
are a set of independent variables. Introduce one hopping parameter for each direc- 
tion, for each plaquette orientation and let G = dn^/ d\ogaQ and H = Bhq/ dlogak- 
Now locking ai = 02 = 03 = a^, the following relations hold at the isotropic point = at. 

dni^/dat = H dKk/dus = G + 2H 

duo/dat = G dKo/dag = 3H. (55) 

Similarly, if 5 = 3/3ofc / d log ao and U = d(3ki / d log ao , then in the ai — a2 — — 
theory, we have at the isotropic point Os = at ^l^ : 

Pok/d log at = S pok/d log as = S + 2U 

Pki/d log at = U Pki/d log as = 2S + U (56) 

In this way the number of functions has been reduced by a factor two. However, to our 
knowledge the derivatives U, S, G, H have not yet been determined individually in the 
A^f = 2 theory. 

By writing the renormalization group equations -r^— = 0, and taking suitable 
linear combinations thereof, one recovers Eq. 2H1 (or more generally Eq. I52j) under the 



dE 

dE 
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consistency conditions 



2{S + U) 



dp 



dlog a 

and obtains the independent sum rule 



k ^ k k<l 



(57) 

3ht-h,|*|[l + 0(a)]. 

(58) 

The vacuum matrix elements vanish in this case by space-time symmetry. 

In thermodynamic studies of QCD using the Wilson regularization |21j, the energy 
density e and the pressure p are extracted in the so-called derivative method by comput- 
ing derivatives with respect to (3 and k of the partition function. Since Z = Tr {T^°}, 
where the temperature T is equal to L^^ , one finds that aV{e — 3p), the 'interaction 
measure', is given by a thermal average (with Boltzmann measure e~^/"^) over the eigen- 
states of the transfer matrix of expression (|52|1. Similarly the expression corresponding 
to aV{e+p) is Eq. (f^ . 



4 Conclusion 

We have extended the known set of lattice sum rules by including the effects of Wilson 
fermions. Our derivation in the transfer matrix formalism kept track of discretization 
errors. The main results are Eq. (|52j) and (|58j) . 

We have also presented new applications of the sum rules. Since they allow us to 
compute the derivative of an energy level with respect to the lattice spacing, we have 
proposed to exploit this to reduce cutoff effects on the spectrum and on the static po- 
tential. One may include the information on the slope in the continuum extrapolation, 
if one disposes of data at several lattice spacings. Alternatively, one can form an im- 
proved estimate of the continuum limit from data at a single lattice spacing under the 
assumption of a functional form for the leading discretization errors, such as the one 
predicted by Symanzik's effective theory |21j . 

Of course these applications are only of practical interest if the three-point functions 
can be computed accurately in numerical simulations. This can probably be achieved for 
the static potential in the pure gauge theory using a suitable multi-level algorithm |25j . 
It also seems realistic in A'^f > 2 theories in the pseudoscalar sector, although high statis- 
tics and efficient all-to-all techniques will be needed. We remark that higher order sum 
rules may be derived straightforwardly by taking additional derivatives of the trans- 
fer matrix with respect to the bare parameters, but the n-point functions are bound 
to become more difficult to evaluate numerically with increasing n. It is nevertheless 
conceptually pleasing that locally the lattice artefacts can be determined at fixed bare 
lattice parameters. 

I am happy to thank Karl Jansen and Rainer Sommer for reading and commenting 
on the manuscript. 
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